BERRY PHASE EFFECTS IN THE DYNAMICS OF DIRAC ELECTRONS 
IN DOUBLY SPECIAL RELATIVITY FRAMEWORK 



Pierre Gosselin 

^ '. Institut Fourier, UMR 5582 CNRS-UJF, UFR de Mathematiques, 

Universite Grenoble I, 
BP74, 38402 Saint Martin d'Heres, Cedex, 
France. 



O 

(N 



O 



o 



X 



Alain Berard, and Herve Mohrbach 



Laboratoire de Physique Moleculaire et des Collisions, ICPMB-FR CNRS 2843, 

Universite Paul Verlaine-Metz, 
57078 Metz Cedex 3, 
France. 

m ' 
> 
ON 

g ! Subir Ghosh 

q ■ Physics and Applied Mathematics Unit, 

C""- . Indian Statistical Institute, 



203 B. T. Road, Calcutta 700108, 
India. 



Abstract: 

We consider the Doubly Special Relativity (DSR) generalization of Dirac equation in an ex- 
ternal potential in the Magueijo-Smolin base. The particles obey a modified energy-momentum 
dispersion relation. The semiclassical diagonalization of the Dirac Hamiltonian reveals the in- 
trinsic Berry phase effects in the particle dynamics. 
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Evidence (see [1] for discussion and references) of ultra-high energy cosmic ray particles 
that violate the Greisen-Zatsepin-Kuzmin bound have compelled theorists to generalize the 
conventional energy-momentum dispersion law of particles, 

p 2 = m 2 , (1) 

based on principles of Special theory of Relativity (SR). 1 The extension requires another ob- 
server independent dimensional parameter, apart from c, the velocity of light. The second 
parameter (k) is expected to be related to Planck energy. Based on this idea Amelino-Camelia 
[2] has pioneered an extended form of SR, popularly known as Doubly Special theory of Rel- 
ativity (DSR). The effect of k appears in the explicit structures of Lorentz transformations in 
DSR, which are non-linear for momenta and momenta-dependent for coordinates. At the same 
time k induces a Non-Commutative (NC) spacetime structure [3], which is referred to as the 
/t-Minkowski spacetime [4]. Generally the scale of k is associated with Planck energy, and all 
K-induced modifications smoothly disappear in the low energy sector (or equivalently in the 
limit k — > oo). In the present paper we take A = (1/k) as the NC parameter and for A = one 
recovers the commutative limit. 

A particular form of DSR extension of (1) that we will be concerned with in the present 
paper is given by, 

p 2 = M 2 [l - \E} 2 , (2) 

with E being the particle energy. This form is known as the Magueijo-Smolin (MS) construction 
[5] of DSR. The corresponding ^-Minkowski NC spacetime [4] is, 

{x\x } = i\x i ; {x\x j } = 0. (3) 

We will comment on the classical nature of NC spacetime algebra, that is being considered here, 
a little later. In the present article we will focus on the DSR generalization of Dirac equation 
[6] in MS base [5]. (For alternative constructions of DSR Dirac equation see [7].) We will take 
up directly the issue of (semi-classical) quantization of the DSR Dirac particle by exploiting the 
generalized Foldy-Wouthuysen (FW) approach [8]. This extension amounts to a semiclassical 
(0(h) quantization of the Dirac particle in the presence of external interactions and in previous 
works of some of us [9] this idea has been successfully employed to reveal intrinsic Berry phase 
[10] effects in the particle dynamics (see [11] for a review). The appearance of Berry phase 
effects is very natural in FW formalism. The Berry potential is induced by the electron spin 
coordinates which are treated as "fast" variables as compared to the position coordinates which 
are considered as "slow" degrees of freedom. The same phenomena was discovered earlier [12] 
when molecular dynamics was studied in the Born-Oppenheimer approximation. In this context 
let us emphasize the novelties in the present work: The DSR Dirac equation will be derived 
in an algebraic formalism which is mathematically very simple and it exploits ideas studied 
in detail in a previous work of one of the present authors [13]. Interestingly this formalism 
allows us to introduce interactions in the DSR scenario in a consistent way. Indeed interactions 
in a DSR framework generally have not appeared in the literature. Lastly the study of DSR 

1 Howevcr, it should be pointed out that very recent data from the Auger cosmic ray observatory does not 
quite support the observation of ultra-high energy cosmic ray particles although results on the ultra-high energy 
cosmic photons is still awaited. 
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Dirac model in FW scheme is entirely new and (not surprisingly) the results obtained by us 
can indicate new directions in this area. 

The paper is organized as follows: In Section II we derive the DSR Dirac equation in 
Magueijo-Smolin (MS) [5] base. Section III is devoted to the FW analysis of the DSR Dirac 
equation constructed in Section II. In Section IV we introduce the external interaction and re- 
veal the intrinsic Berry phase contribution in the present case. The paper ends with conclusions 
and areas of future studies in Section V. 



DSR DIRAC EQUATION IN MAGUEIJO-SMOLIN BASE 



Let us first put our approach in its proper perspective in the context of quantum NC DSR 
theories. Our aim is to exploit the canonical framework [13] in constructing the theory and 
subsequently consider its quantization. For this reason our analysis is completely classical, 
at least for the time being. Thus we use (classical) Poisson brackets in place of (quantum) 
commutators. On the other hand, the authors of [14] discuss a quantum DSR theory from the 
very beginning and consider the spacetime as quantum in nature. However, the classical nature 
of the field variables are retained in [14] in the sense that products of fields are not replaced 
by their *-product. In this section we will construct the DSR Dirac equation in MS base (see 
Appendix). To motivate the present derivation, we have to introduce the full ^-Minkowski 
Non- Commutative (NC) phase space, 

{x\x } = Xx l ; {x\x j } = ; {x\p>} = ; {p»,p»} = 0. (4) 

{x ,^} = Xp l ; {x\p } = ; {x°,p } = -1 + Xp°. (5) 

We are in the classical framework and will interpret the phase space algebra as Poisson brackets. 
Our metric is diag g 00 = —g u = 1. For convenience it is expressed in a covariant form, 

{x^, X v ^ X{x^Tj u X u Tj^j, 

{x^Pu} = -sv + {p^Pu} = o, (6) 

where 770 = 1, 77^ = 0. For completeness we mention that the Lorentz generator 

satisfies the undeformed Lorentz algebra, 

{J"", J a(3 } = g^J va + g m J pv + g v Pj aix + g va J^, (8) 

but induces deformations in the transformation laws. The energy-momentum dispersion law, 
consistent with the form of given in (7) is the MS relation (2), 

p 2 = M 2 [l-XE] 2 = M 2 [l-X(r ] p)} 2 . (9) 

Indeed in all the above computations one uses the NC phase space algebra (6). 

It is quite obvious that a direct generalization of the physical laws in the NC phase space is 
complicated, both mathematically as well as conceptually [6, 7]. However, there is an easy way 
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out [13] which we now explain. We can introduce a map (X^P^) = F~ 1 (x f j,,p fJ ) [13] between 
(xpjPn) - the physical /t-Minkowski NC phase space and (X M , P M ) - a completely canonicalpha.se 
space. Explicitly the transformation rules are the following: 

The (X M , P M ) variables obey canonical Poisson Bracket algebra 

{X li ,P v } = -g lu/ ; {X fl ,X u } = {P fl ,P u } = 0. (11) 

We have also shown in [13] that (X^, P v ) have conventional (Special Theoretic) Lorentz trans- 
formation properties whereas, as mentioned before, {x^p v ) obey /t-deformed Lorentz transfor- 
mation laws [15]. 

The above relations in (10) are invertible, 

,„ = X„(l + A(„P)) = + AP ); P „ = ITT ^ y = JT ^ p - y (12) 

We note that the above relations are classical and one has to order them appropriately under 
quantization. However, we will show these mappings survive under some restrictions and in 
the present work that is sufficient. 

Our framework of generalizing physical laws to /t-Minkowski phase space is the following: 
Start with the known form of a relation in the (auxiliary) canonical phase space (X^, P u ). Now 
simply map the relation using (10) to the physical NC A-variables (x^,p u ). This yields the 
cherished A-deformed physical law. Indeed it is not obvious that the procedure will work but 
we have explicitly shown its validity in various instances in [13] in the context of point particle 
models. We will demonstrate (see Appendix) that this principle works for the Dirac equation 
as well. 

To this end, we start with the momentum space Dirac equation in canonical phase space: 

(7 M -P/x - M)u(P) = 0, (13) 



{l" - \ p) , - M)u{p) = 0. (14) 



where 7^ are the standard "7 matrices". We now provide the one line derivation of the DSR 
Dirac equation in MS base: 

Pjj_ 

(I — (ZZE; 

V K 

This is the A-extended Dirac equation in MS basis. This is a new result. In the Appendix, we 
substantiate the validity of our claim by comparing (14) with existing results in the literature 
[6]. 

For applying the FW transformation on our MS-Dirac equation (14) we express it in Hamil- 
tonian form, 

(E-a.p- (3M(l- \E))u = 0. (15) 
After rearranging (15) we obtain, 

E(l + XM(3)u = (a.p + (3M)u (16) 

Eu = { \Z^ (S.p + (3M)n. (17) 
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Quite surprisingly we find that the Hamiltonian H (or E in (17)) has a non-hermitian structure: 

H = A + B(3 + Cap, (18) 

where 

-AM 2 



A = 
B = 
C = 



1 - A 2 M 2 
M 

1 - A 2 M 2 
(1 - AM/3) 
1 - A 2 M 2 ' 



(19) 
(20) 
(21) 



This is a new result. This feature was not revealed in previous analysis [6, 7] since a proper 
quantum treatment of the DSR Dirac equation was not attempted. This problem is tackled by 
performing a similarity transformation on H with the matrix D, 



D = C- 1 ' 2 



1 + —p = a + b{3, 

K 



where 



° = ^2 



M 2 



K 



This leads to the final Hermitian form of H: 



M 2 



H = D(A + BI3 + CapjD- 1 = A + Bj3 + (1 - ^-)-*a.p 



(22) 



(23) 



(24) 



This H is suitable for the conventional FW procedure with similarity transformation performed 
by the unitary matrix U , 



U 



1-A 2 M 2 



+ B' 



1~X 2 M 2 



+ B 2 + B 



1/2 • 



(25) 



which can also be written as, 

(^(1-A 2 M 2 )^ + M 2 + M + VI - \ 2 M 2 (3a.p) 



U 



(2tJ (1 — A 2 M 2 ) p 2 + M 2 (yj(l- \ 2 M 2 )p* + M 2 + M 



1/2 



(26) 



In this way we see that U is very similar to the usual FW transformation U for a free particle 
of momentum ~j? 2 which is 



U 



{^J^ 2 + M 2 + M + fia.p^j 
2^ 2 + M 2 {^J^ 2 + M 2 + m) 



1/2- 



(27) 
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Indeed, U in (26) reduces to U for A = 0. 2 The energy eigenvalues of H are: 



/t 2 1-A 2 M 2 VU-A 2 M 2 / V1-A 2 M 2 



, A2 |/2 , AM 2 ± V(l - A 2 M2) f 2 + M^j (28) 

Clearly they satisfy the MS dispersion law (2). This can be considered as an alternative 
way of obtaining the MS dispersion law starting from the MS extension of Dirac equation, in 
the spirit of [18]. We emphasize the framework adopted by us, i.e. constructing the MS Dirac 
equation by exploiting the canonical phase space coordinates and subsequently utilizing the 
FW formalism to compute the energy eigenvalues that are identical with the MS energy values, 
is quantum mechanical in nature and is totally new. 

We would like to make a comment on the particle and anti-particle spectra corresponding 
to the MS dispersion law, although this is not directly relevant for the present work. For 
the normal Dirac particle in commutative spacetime, the particle and anti-particle spectra, 
E± = ±y/p 2 + M 2 are symmetrically placed with respect to the zero-energy level. This is 
obtained by putting A = in (28). At the same time it is evident that this feature is absent 
for the MS spectra (28) with a non-zero A. However, we can restore this symmetry simply by 
shifting the zero level of the energy value by the constant amount 1 Z\2 4 M 2 in (28). Thus the 
particle and anti-particle sector energy levels are given by, 

— X 2 M 2 )fP + M 2 

E± = ± — — • (29) 



PHYSICAL POSITION OPERATOR AND BERRY CURVATURE EFFECTS ON 

PARTICLE DYNAMICS 

It is well known that, (in the normal commutative spacetime), the position operator r in 
Dirac equation does not represent the physical coordinate of the particle simply because the 
magnitude of the particle velocity obtained from taking time evolution of r is c, even for a 
massive particle. On the other hand the FW formalism provides a natural prescription for 
a physical position operator r + that correctly reproduces the particle velocity. Incidentally 
r + operators are noncommutative in nature and as we shall see below, this feature can be 
attributed to a Berry phase effect. Phenomena of a similar have been observed before [19] in 
condensed matter systems. In the present case we find that the Berry curvature effect is further 
modified by the A-corrections. 

We now limit ourself to the dynamical evolution of a positive energy particle whose energy is 
formally the projection (we denote the formal projection operation by V) of the Hamiltonian on 
the positive energy subspace. Physically, we thus consider the adiabatic approximation which 
allows to neglect the interband transition, by identifying the momentum degree of freedom as 
slow and the spin degree of freedom as fast, similarly to the nuclear configuration in adiabatic 

2 Once again, in an algebraic way, U can be computed directly from Uo by the change of variable ~]t — > 
\f\ — \ 2 M 2 p. This is simply because in our equation (24) the factor (1 — A 2 M 2 )~2 appears in front of a.p. 
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treatment of molecular problems. To be coherent with this projection on the positive energy 
subspace, the same projection has to be done for all dynamical operators. Therefore the physical 
position operator ~r* + is obtained as 

1> + =V (UD^D- 1 ^) = V (UT>U + ) (30) 

with the MS position operator r — %h^. Notice that we are using the canonical representation 
of the position operator, (valid for commutative spacetime), for the A- variables as well. This 
is allowed if one recalls from (4,5,6) that the spatial sector remains unaffected in this NC 
spacetime. This can also be justified in a more rigorous way by going to the canonical phase 
space framework (10). This leads to the following explicit form of r+: 

( d \ d ( dU + \ 

r + =i%v u^—u + =ih— + v uih—^ \=r + A, r) (31) 

\ op ) op \ op J 

where we have introduced Am, the so called Berry connection. The (r) in Am indicates that 
this Berry connection accompanies the position ~r\ (In general, Berry connections can appear 
both in coordinate and momentum.) Indeed the state space of the Dirac electron is spanned 
by the basis of plane waves of the form \p,<p a ) = \p) \ fa), where \ip a ) is a zero-momentum 
spinor, such that U(p) \ip a ) = \u a (p)) which is a spinor solution of the Dirac equation. For 
\<f a ) = (1, 0, 0, 0) and |</?g) = (1, 0, 0, 0) the matrix elements of Am are then 

(Ar)) =ih(<fa\U (p) \tp p ) = ih (u a (p)\ -—^ \up(p)) . 

v 7 «P op op 

This is the definition of the Berry connection. The explicit structure of Am in the present 
context is, 

A( r ) =V[U rU 1 = h — v -, (32) 

V ' 2^(1 - \ 2 M 2 ) f + M 2 (^(1 - A 2 M 2 ) p 2 + M 2 + M J 

with "o^ being the Pauli matrices. Note that for a free particle the momentum is FW invariant, 
i.e. J? + =V (UjfU^ 1 ) = ~p*. 

As we have mentioned before, the physical position operators r + are no longer commutative: 

r\y + ] = ihe^Q k = l h((v\^A\ r) -V\^ 
[p\jP] = 

p\ri} = -ihS ij (33) 



with the non-abelian Berry curvature defined by 

3 = h - Jm* + , (l-^)OT V (34 ) 

((1 - \ 2 M 2 ) p 2 + M 2 f 2 \ - \ 2 M 2 ) p + M 2 + MJ 

The dynamics is generated by the Hamiltonian equations of motion, 

= i[H,0], (35) 
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where O is a generic operator and H is the Hamiltonian. In the present instance with E + 
representing H, we find, 



I = = ; ^=0. (36) 

y/(l - \ 2 M 2 ) -f 2 + M 2 

Also one finds that the non-relativistic expression for velocity is obtained for M — 1/A which in 
fact corresponds to the upper bound of mass (or energy) available for the particle. In the next 
section we will study the dynamics of the MS Dirac particle in presence of external interaction. 



DSR DYNAMICS IN PRESENCE OF SCALAR INTERACTION AND BERRY 

CURVATURE EFFECTS 

In general it is not possible to perform the FW transformation exactly when interactions 
are present but one can always resort to a perturbative study for a " small" potential term. But 
here instead of a perturbative expansion we consider a semiclassical approximation (0(h)) and 
in addition limit ourself to 0(A) effect only. 

Once again we start with the free MS Dirac equation (14) in canonical phase space and 
introduce an interaction potential V(R), 

{^Pfx - M - 7 V(R)}u(P) — 0. (37) 

This interaction can be thought of as scalar component of U(l) interaction, in a particular 
frame where the vector potential vanishes. 

Translated to A-variables this equation reads, 

- P \( V p)) ~ M ~ 7V(f(1 _ Ai ^H£) = 0. (38) 

Let us consider the A in its lowest non-trivial order. Hence to 0(A) the Dirac equation is, 

(1 + XV (T>) +A r\W (-t)) E u(p) = (A + B'(3 + C(a.p) + V (t))u(p). (39) 

with B' = B- XMV ( r>) and A, B, C given in Eqs.(19-21) 
To fit in the FW scheme (39) is rewritten as, 

Eu(p) = (1 - XV ( -t) - X r>.W ( -t)) (A + B'(3 + C(a.p) + V (^))u(p). (40) 

The same procedure as in the free case is needed to render this H hermitian and we obtain a 
Hamiltonian that is hermitian in O(A) and h: 

H = D(1-XV (r^-A^.VU (!>))(A + B'f3 + Ca.p + V (^/T 1 
= (1 - XV (T) - XT.VV (!>)) D(A + B'(3 + Cap + V (^D- 1 

= {1 - W (-?) - X? .W (-?)) Ia + B'(3 + (1-— )-*a.p. + V(-t)\ (41) 
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In ref [20] we have shown how to diagonalize Dirac like Hamiltonians in the presence of external 
fields at the semiclassical order ( to get corrections beyond the semiclassical order see ref 
[21]). During this process of diagonalization, we have shown that in general (in the presence 
of electromagnetic or gravitational fields), both position and momentum operators acquire 
a Berry-phase contribution making the coordinate and momentum algebra noncommutative. 
Note that in the presence of interactions, the projected (on the positive energy subspace) 
dynamical operators emerge naturally during the diagonalization, without resorting to the 
adiabatic approximation. 

Therefore applying this diagonalization procedure in our case leads for the positive energy 
sector to the following semiclassical diagonal Hamiltonian 

H + = V (UHU- 1 ) (42) 

where U is the same matrix as in the free case with M replaced by an effective mass M (1 — XV ("f*)). 
This leads to 

H + = (1 - XV ( r>+) —\T* + W ( T+)) (^/t 2 + (l-XV ( r>+)) 2 M 2 - AM 2 + V ( r>+)) 

+o(\J\ 2 + h 2 ) (43) 

where ~r* + = "r* + «4(A,r) with «4(A,r) the same as in the free case with the effective mass 
M (1 — XV ("f > + )) . We observe a renormalization of the energy through the term 

(1 - (T>+) -A rVW ( • 

Note that the momentum get also a Berry potential correction which at the order considered 
is negligeable = J? + 0(Xh). The commutation relations are also unchanged exept for the 
mass renormalization. Hence, the dynamics in presence of the potential is derived to be, 

t + = V(p)H + - f X0 ; ^=-V (r+) if+, (44) 

where is given in (34) (with M replaced by an effective mass M (1 — XV (V')) ). Note that 
the effect of Berry curvature is manifested very directly once an interaction is present. Indeed, 

the equation for the velocity contains an anomalous velocity term ~]f x© of order H which 
causes an additional displacement of the electrons orthogonally to the momentum J? . This 
phenomenon depends on the particle spin through the Berry curvature. Typically, in canonical 
spacetime, with V being an electrostatic potential, this is the well studied spin Hall effect. In 
fact it should be mentioned that the topics of spin Hall effect of electrons or photons, gravita- 
tional Hall effect for photons, that have become areas of intense research, all owe there existence 
to this type of Berry curvature effect. However, as we point out below, modelling the analogue 
of spin Hall effect for MS particle is probably more complicated. 

CONCLUDING REMARKS 

In this paper we have considered particle dynamics in /t-Minkowski spacetime which is 
a particular form of noncommutative spacetime. Indeed quantization of the particle model 
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is tricky because of the operatorial form of noncommutativity involved in the phase space 
commutation relations. The novelty of our work lies in the fact that we have been able to 
bypass this problem by exploiting a semiclassical approach of the Foldy-Wouthuysen formalism. 
Indeed this makes our analysis semi-classical in nature with quantum effects in {0(h) taken 
into account. We then have shown that the dynamics of MS particles in the presence of an 
external field should be influenced by a the ^-parameter but also by the presence of a Berry 
curvature which is itself k dependent. 

In this connection, let us pause to mention the significant (and possibly more ambitious) 
work [22] that attempts to deal with the full ^-Minkowski quantum field theory. Considering 
the *-product for the fields, the authors of [22] demonstrate that if properly interpreted, the 
interaction vertices have the requisite symmetry under the interchange of momenta of identical 
incoming particles. Indeed, the ^-Minkowski *-product is much more involved that (Moyal) 
♦-product that appears when the non-commutativity is constant in nature. It is true that 
we have so far only looked at the quantum mechanics of a DSR model but the nature of the 
energy-momentum conservation laws [13] inherently possess the above mentioned symmetry. 
This seems to indicate that the DSR in Magueijo-Smolin base might be a better alternative 
that DSR in other bases in the context of formulating the DSR quantum field theory. 

The next task we wish to pursue is the effect of full electromagnetic interactions on the DSR 
particle model considered here. This is a non-trivial extension since electromagnetism involves 
a U(l) gauge symmetry and one has to appropriately generalize the concept of gauge invariance 
in noncommutative spacetime. This feature probably did not show up in the present restricted 
setup. 

APPENDIX: 



We will compare the MS Dirac equation obtained here in (14) with that derived in [6]. It 
is to be noted that [6] uses the bicrossproduct base whereas we have used the MS base. 

It is well-known in the DSR community that there are distinct formulations (or bases) of 
the DSR regarding the modified dispersion relation, transformation properties and phase space 
algebra. The most popular are known as bicrossproduct base [4] , MS base [5] and standard base 
[16]. The different bases are connected by non-linear transformation and [17] but the different 
bases are inequivalent with drastically different physical consequences due to the non-linearity 
in the transformations. 

The DSR Dirac equation obtained in [6] is 

( 7 ^ - u)u = 0, (45) 

where, 

e XE ' - cosh(Acu) , . 

a = . e " m . (47) 

k smli (Aw J 

Here (£",pQ are energy and momenta in bicrossproduct basis and they satisfy the dispersion 
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law [4], 

-^cosh(AE') - (p) 2 e XE> = 2K 2 cosh(\u). (48) 
A 2 

The relations connecting the above set to the MS set of variables (that we have used) are [17], 

p'. = Pl - E' = -\log(l - 2AE + X 2 f),E = -\-(l - e - 2XE ' + \ 2 fi 2 ) ; Am = tanh(Xuj). (49) 
2 A 2 A 

It is straightforward to check that using (49) on (45) will generate (14), the MS Dirac equation 
derived in our paper. 
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